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THE MUON CONTENT OF EAS AS A FUNCTION OF PRIMARY ENERGY 

P. R. BLAKE, W. F. NASH, M. S. SAICH and A. J. SEPHTON 
University of Nottingham, England. 


1 ' Introduction At Haverah Park the muon content of EAS has been 
measured over the wide primary energy range 10 1 to lCT^ eV. At the 
Bangalore Conference it was reported (Blake et al, 1983) that the 
relative muon contey^ of E^S decreases smoothly [~ P(600 m) ’® 4 ) over 
the energy range 10 - 10 eV and therefore concluded that the 

primary cosmic ray flux has a constant mass composition over this 
range. In the same paper (on the basis of a preliminary analysis) 
it was also reported that an apparent significant change in the 
power index occurs below 10 eV [ p (250 m) * U ] . Such a change 

would indicate a significant change C in primary mass composition in 
this range. 

Further anaj.^sis confirms the earlier conclusions concerning EAS of 
en ?§^y - 19 eV * However post- Bangalore analysis of data in the 
10 - 10 eV range revealed a previously overlooked selection bias 

in the Bangalore data set. This paper presents the full analysis of 
the complete data set in the energy range 10 - 10* eV with the 

selection bias eliminated. 

2. Methods Three 10 m 2 muon detectors (-* 400 MeV threshold) were 
operated in conjunction with the deep water Cerenkov tank array 
(operated by the University of Leeds EAS group) at the Haverah Park 
site. Two of these muon detectors were immediately adjacent to 
large area water Cerenkov detectors, enabling a direct measurement 
of the ratio of the muon response to the water Cerenkov response 
(P M /P c ) t0 be made - Th ® water Cerenkov response ( P ) is dominated 
at low core distances by the electromagnetic response. The directly 
measured ratio (P/P ) has the advantage of not being very sensitive 
to EAS parameters and therefore not sensitive to any inaccuracies in 
determining these parameters. 

, 3 * Results (10 - 10 17 eV) The basic problem in establishing an 

accurate determination of the muon content over a large primary 
energy range arises from the measurement technique. Inevitably in 
most arrays (including the Haverah Park array) the core distance 
range covered is a strong function of primary energy (or shower 
size). The ratio (P/P c ) itself will be a function of core distance 
and thus direct determination of p /p at a specific core distance 
is limited to a small ^range of primary energies. If a relationship 
of the form p = k is adopted, k can normally be regarded as 

constant at a H particular core distance (R) ; the power index 
* (normally 0.94) may be a function of core distance and shower size. 

A core distance dependence would arise if the lateral density 
distributions of p or p significantly change with shower size. 

„ Measurements indicate that the lateral distribution of p is 

unchanging at energies > 10 eV. However the lateral diitribution 
function of the electromagnetic component (and therefore p ) does 
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show some change as the electromagnetic cascade maximum penetrates 
more deeply into the atmosphere as E increases. Thus in 
determining the exponent ' a * it is slnsible to limit the core 
distance range to as small a core distance interval as is 
statistically feasible. 


The available experimental data were divided into three overlapping 
intervals of core distance:- 50 m <_ R < 150 m? 100 m <_ R <_ 200 m? 
and 150 m £ R _< 250 m. The data were further subdivided into four 
sec 0 intervals 1 -* 1.1 -* 1.2 "*1.3 ”*1.4 and six shower ' log ' size 
intervals. The data were then used to calculate the mean measured 
muon and water Cerenkov responses at R = 100 m, 150 m and 200 m. The 
results for 1.0 < sec 0 < 1.1 are displayed in Table 1 & Figure 1. 

Clearly the p /p 
ratios change^ C 
smoothly over 
this decade of 
primary energy. 

Table 2 shows the 
best fit values 
of a obtained 
from linear 
regression fits 
on the 

relationship p = 
kp a for all cHe 
R and sec0 
intervals . 



FIGURE 1 
distances 


Measured p versus measured p at three specific core 
M c 


R 


/ m 2 
P 1 50 

N 

/m 2 

% 

+ Ap 
“ H 

/ m 2 

vv 

+ Ap /p 
— |i c 

100 

m 

1.23 

2466 

0.560 

0.011 

3.306 

0.169 

0.003 



1.81 

2362 

0.805 

0.017 

4.833 

0 . 167 

0.004 



2.66 

1793 

1.188 

0.028 

7.091 

0.168 

0.004 



3.90 

1106 

1.664 

0.030 

10.586 

0 . 157 

0.003 



5.73 

529 

2.430 

0.11 

15.841 

0.154 

0.007 



8.41 

234 

3.764 

0.25 

24.209 

0 . 155 

0.010 

150 

m 

1.23 

2040 

0.306 

0.007 

1.333 

0.230 

0.005 



1.81 

2303 

0.432 

0.009 

1.915 

0.226 

0.005 



2.66 

2065 

0.634 

0.014 

2.773 

0.229 

0.005 



3.90 

1397 

0.930 

0.025 

4.119 

0.226 

0.006 



5.73 

820 

1.301 

0.045 

6.239 

0.209 

0.007 



8.41 

378 

2.025 

0.104 

9.907 

0.205 

0.010 

200 

m 

1.23 

904 

0.179 

0.006 

0.648 

0.276 

0.009 



1.81 

1335 

0.272 

0.007 

0.979 

0.278 

0.007 



2 . 6,6 

1540 

0.368 

0.010 

1.438 

0.268 

0.007 



3.90 

1287 

0.584 

0.016 

2.098 

0.278 

0.008 



5.73 

886 

0.809 

0.027 

3.255 

0.249 

0.008 



8.41 

470 

1.255 

0.058 

5.141 

0.244 

0.011 
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Table 1 (previous page) Mean muon (~) and water Cerenkov pH 

densities as a function of cfere distance (R) and shower 
size (p g) - 1.0 ^ sec0 £ 1.1. 

[ P 1 50 water Cerenkov response in vertical equivalent 

muons at 150 m from the EAS axis. N = no. of EAS used in 
the analysis] 


R 

Sec0 

1.0 - 1.1 

See© 

1.1 - 1.2 

Sec0 

1.2 - 1.3 

Sec© 

1.3 - 1.4 

100 m 

0.95 (0.01) 

0.92 (0.02) 

0.90(0.02) 

1.01 (0.04) 

150 m 

0.94 (0.02) 

0.95 (0.02) 

0.99(0.02) 

0.99(0.06) 

200 m 

0.94(0.02) 

0.94(0.03) 

0.93 (0.08) 

0.90(0.05) 


Table 2 Values of a in the relationship p 
R and ©. 



as a function of 


Clearly, at least for low zenith angles, an average value of a = 
0.94 fits the data well in the core di^ance range 50 m m and 

over the primary energy range 1.6 x 10 eV < E < 1.6 x 10 eV. 


4 . Muon^onten^ lO 


16 


19 

- 10 eV 


Data acquired in the primary energy 
range 10"' - 10* w eV from the infilled 500 m detector array at 
Haverah Park enables p (200 and P (200 m) to be determined for 

the decade of energy afeove 10 eV. ¥he resul^g are jotted in 

Figure 2 which shows the complete data from 10 - 10 eV. 


E P [eVl 

io’ 6 10 17 10 M 



FIGURE 2 Measured p versyg p^ a^g200 m over energy range 
approximately 10 - C 10 eV. 
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It is clear that a relationship of the form p (200 m)ocP c (200 m) 
fits smoothly and well over these two decades of energy* This 
implies a smooth gradual change in the height of the electromagnetic 
cascade maximum throughout these two decades. 


Thus the muon data yields strong eviden^g that |.^ere is no marked 
primary mass change occurring in the 10 to 10 eV region as has 
been suggested. If there is a significant change in primary mass, as 
suggested by Chantler et al (1^3), then the muon density would be 
approximately 50% higher at 1<J eV than is measured leading to a 
value of a ~ 0.8 [Blake et al, 1983]. Such a low value of a is 
clearly incompatible with the measured data (see Figure 2). 

17 18 1 8 
Comparison of the 10 to 10 eV decade with energies above 10 eV 

has been carried out at a core distance of R = 600 m. In this 

region the data and conclusions remain unchanged since Blake et al 

(1983). Again the observations yield a^gmooth and normal change in 

muon content in the range 10 eV to 10 eV, implying no 

significant primary mass change over the energy range. 

5. Discussion and Conclusions The results presented above lead to 
the conclusion that the muon content of EAS (relative to the water 
Cerenkov response) decreases smoothly in a constant manner at 
primary energies above 10 eV and therefore suggest a constant 
primary mass composition in this energy region. A direct and 
constant relationship between the primary mass energy .(E ) and the 
water Cerenkov response p (R) is assumed. The water Cerenkov 
response at 500 m core distance, p (500) , [or at 600 m,p(600)] is 
used as the shower size parameter at Haverah Park for energies > 

10 1 eV and the relationship between p (500) and E is of thg ^rm E 
= 3.87 x 10 Ipl500)] [HILLAS, Model A]. Thfis p p Q 

implies p ocE at 500 m from the EAS axis. 

M P 

Assuming therefore that this result does indj.gate a constant primary 
mass composition at energies greater than 10 eV it does give 
strong support to the results from other studies of muons in EAS at 
these energies. The Akeno group (Nagano et al, 1984) find that t^g 
size spectra for both electrons and muons are unchanging above 10 
eV. Thus any change from heavy primary nuclei ty^a light mass 
composition flux must occur at energies below 10 eV. 
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